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A  METHOD  TO  ESTIMATE  THE  VERTICAL  DISPERSION  PARAMETER  IN  A  10  Km 
RANGE 

Lei  Xiaoen,  Jia  X inyuan,  and  Yang  Jinte 
(Institute  of  Atmospheric  Physics,  Academia  Sinica) 

ABSTRACT 

Based  on  the  Monin-Batchelor  Similarity  Theory  and  the  concept 
of  effective  roughness  length,  this  paper  presented  an  empirical 
vertical  dispersion  model  in  a  10  kilometer  range.  It  could  be 
used  under  a  flat  and  homogeneous,  as  well  as  complex, 
topographical  condition. 

I .  FOREWORD 

The  application  of  the  similarity  theory ^  to  a  flat  and 

uniform  terrain  in  a  close  distance  (within  2  kilometers)  has 

(2  3 ) 

already  had  precise  conclusions'  *  .  But,  what  happens  when 

it  is  applied  to  the  situation  of  a  complex  terrain  and  a  further 

distance  (for  example,  as  far  as  10  kilometers)?  This  is  a 

problem  of  great  concern  to  us. 

In  this  paper,  the  concept  of  similarity  theory  was  applied 

to  derive  the  vertical  dispersion  model  (i.e.  the  correlation 

between  the  vertical  dispersion,  parameter  £,  and  the  downwind 

distance,  X)  for  a  flat  and  uniform  terrain.  Next,  the  results 

(4) 

were  compared  to  the  turbulent  dispersion  model  which  is 
widely  used  in  the  plain  area.  Its  applicability  in  a  10 
kilometer  range  was  discussed.  Thirdly,  the  concept  of  effective 
roughness  length ^  was  used  to  extend  the  dispersion  model  from 
a  plain  terrain  to  a  more  complex  topography.  Furthermore, 
through  a  comparison  with  the  observed  data,  the  effectiveness 
and  applicability  of  this  dispersion  model  were  discussed  under 
a  complex  tower  surface  condition. 

II.  DERIVATION  OF  THE  CORRELATION  BETWEEN  Z  AND  X 
According  to  the  method  of  Ito^,  and  usual  expression 

between  5?  and  z  is 

*-  £js(i)f(C)C'*  +  A , 
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where  u(z)  is  the  average  wind  speed ,  u*  is  the  speed  unit  on 
the  ground  level,  Z  *  %,  L  is  the  length  unit  on  the  ground 
level,  f ' (?)  *  f(c)  is  a  generalized  function,  A  is  a 

constant  to  be  determined  and  k  is  the  Kalman  constant.  Both 
were  chosen  to  be  0.36. 

In  the  neutral  situation,  the  logrithmic  wind  profile 
expression  was  substituted  into  (1) .  After  integration  and 
calculation,  one  gets 

*~pK  +  t-]*  <2, 


when  zQ  is  the  roughness  length  of  the  ground. 

In  non-adabatic  situations,  we  substituted  the  logrithmic- 
ally  linear,  BWIB'  type,  and  exponential'  wind  profiles  into 
(1) ,  respectively.  Through  integration  and  rearrangement,  we 
finally  obtained  the  following,  respectively 


*-  -?[£<*' -a+ <*—»(£  -£.*) 
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where  6  is  the  generalized  constant  in  the  logrithmically  linear 

(9) 

wind  profile.  In  later  calculation,  it  was  chosen  to  be  6 

y  ■  (1-r  r)1//4*  Y  *  (1-r  and  r  was  chosen  to  be  15  in 

Ii  O  “ 

computations v  1  .  m  is  a  parameter  related  to  the  stability  and 
topographical  roughness,  and  B  *  7/z0. 


III.  VERTICAL  DISPERSION  MODEL  OF  FLAT  AND  UNIFORM  TOPOGRAPHIES 
In  order  to  derive  the  vertical  dispersion  model,  the  key 
is  to  properly  choose  the  stability  parameters  L  and  m.  The 


2 


so-called  flat  and  uniform  terrain  is  an  idealized  condition, 

which  for  z  is  .01  meter.  In  the  following  calculation,  the 
( 4) 

Pasquill  stability  classification  method  was  adopted.  The 
stability  parameters  L  and  m  were  chosen  based  on  existing 
measured  results^4'^®'^'^^  .  The  values  are  given  in  Table  1. 

Table  1.  Stability  Constants 

1.  Stability  Classification 

2.  L (meter) 


•i  •»» 

:  a  ■ 

*  i 

C 

!  o 

I  *  | 

7 

•  ! 

|  -0.19 

-0.13  ! 

-0.1 

I  -0.02 

- 

- 

J.  U*) 

— 

-  ! 

- 

~ 

130 

33 

After  the  parameters  were  selected,  calculations  were 
carried  out  with  respect  to  (2)  -  (5) .  A  comparison  of  the 
results  showed  that  the  exponential  pattern  described  by  equation 
(5)  was  most  suited  for  types  D-A.  For  types  E  and  F,  it  was 
better  to  use  the  linear  logarithmic  relation  (3) .  Hence,  we 
finally  obtained  the  following  vertical  dispersion  model  for  a 
flat  and  uniform  topography. 

I.  A  35:  X  -  0.15  +  2.09  X  10‘**«  -  1.13  X  10'i**1 

4.  B  £:  X  —  0.13  +  3.02  X  10'***  —  1.84  X  lO'J'" 

3.  C3S:  X  -  0.11  4-  5J9  X  10'i*"  -  3.83  X  lO'i"* 

4.  D 2ft:  X  -  0.09  4-  3.59  X  10'i*"  -  3J3  X  10V*  (6) 

j-,  E&:  X  -  7.69tInJ  4-  2.77  X  10‘#  +  7.S5  X  10”'**  + 

1.54  X  lO-'i’lnl  4-  4.1  X  10“V 
4,  F  5t;  X  -  7.69flai  4-  2.77  X  10»*  4-  3  J7f»  4- 
$.59  X  lIT’lla*  4-  7.54  X  lO-’f* 


1.  Type  A;  2.  Type  B;  3.  Type  C;  4.  Type  D;  5.  Type  E?  6.  Type  F. 

Can  the  dispersion  model  (6)  be  applied  to  a  10  kilometer 
range?  We  carried  out  a  comparative  analysis  between  the  calculated 
results  obtained  from  equation  (6)  and  the  widely  recognized 
Briggs*4*  intrapolation  formula  (applicable  within  10  kilometers)  /370 


3 


in  a  0.1  -  10  kilometer  range.  Within  this  range,  a  total 
of  20  points  had  been  chosen  for  statistical  analysis.  The 
a z  values  in  each  selected  distance  ( az  is  the  standard 
deviation  of  the  smoke  mass  particle  distribution  in  the  vertical 
direction.  Under  a  normal  distribution  situation,  the  relation 
between  z  and  z  is  <?_  *  1.25  z.  This  relation  was  used  in  the 

.  z  (13) 

conversion  in  this  paper'  .)  was  compared  with  the  Briggs 
values.  The  ratios  are  given  in  Table  2. 

Table  2.  Comparison  of  the  Dispersion  Models  in  Eq.  (6)  and 
the  Briggs  Model  in  the  0.1  -  10  Kilometer  Range. 

1.  Stability  Type;  2.  az  (this  work)/az  (Briggs) 


From  the  table,  it  can  be  clearly  seen  that,  within  the  0.1  -  10 
kil^neter  range,  the  dispersion  model  (6)  and  the  Briggs 
intrapolation  formula  are  very  consistent. 

When  the  vertical  dispersion  is  estimated  empirically  by 
using  the  similarity  theory  concept  to  a  range  as  far  as  10 
kilometers,  the  problem  is  whether  a  puff  of  smoke  is  still 

located  in  a  constant  flux  layer.  With  regard  to  the  thickness 

( 8 ) 

of  a  constant  flax  layer,  Monin  pointed  out  that  the  average 

(14) 

is  50  meters.  Thuillier'  also  pointed  out  that  within  a  50 
meter  thick  ground  layer,  the  similarity  theory  determing  the 
mean  wind  speed  profile  is  a  useful  tool.  Hence,  we  can  consider 
that  a  constant  flux  layer  is  at  least  50  meters  thick.  But, 
the  thickness  of  this  layer  has  an  apparent  relation  with 
stability.  Under  a  neutral  condition,  it  is  usually  chosen  to 
be  100  meters,  which  has  already  been  proven  by  a  lot  of 
observation  data*i5^.  Furthermore,  with  increasing  roughness, 
its  thickness  will  further  increase.  As  Soma  pointed  out, 
for  a  city  area,  the  thickness  was  at  least  200  meters.  Lappe 
also  pointed  out  that  its  thickness  can  reach  140  meters. 
Panofsky^1®*  pointed  out  that  under  a  gust  condition,  the  variation 
of  u*  was  still  considerably  small  up  to  an  altitude  of  150  meters. 
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For  an  unstable  atmosphere,  the  vertical  variation  of  u*  is  even 

slower.  Based  on  the  meteorogical  observations  at  several 

positions  along  the  Japanese  coast,  Sensku^7^  discovered  that 

the  logarithmically  linear  region  between  wind  and  temperature 

frequently  exceeded  100  meters.  Sometimes,  it  is  as  high  as  200 
(18 ) 

meters.  We  had  also  discovered  that  the  logarithmically 
linear  law  was  well  satisfied  up  to  an  altitude  of  100  meters. 
Therefore,  for  a  neutral  condition,  the  constant  flux  layer  is 
chosen  to  be  over  100  meters.  Even  up  to  200  meters  in  altitude, 
the  error  would  not  be  too  large. 

In  order  to  explain  the  average  altitude  of  the  smog  when 
the  horizontal  range  is  10  kilometers,  we  presented  the  values 
of  z  and  ^  at  a  10  kilometer  distance  in  Table  3. 

Table  3.  Average  Vertical  Displacement  of  Smog  Particles  and  Z/L. 

1.  Stability  Type,  2.  z(m)  at  a  distance  10  kilometer 


1  un 

A 

1 

C 

D 

E 

F 

lO*****) 

2390 

1230 

MS 

190 

<1 

31 

I/L 

— 

— 

- 

- 

0.41 

O.tt 

From  the  table,  one  can  see  that  the  smog  can  still  be 
considered  to  disperse  in  a  constant  flux  layer  (50  meters  in 
thickness)  within  a  10  kilometer  range  in  approximation  for  types 
E  and  F.  Moreover,  the  values  of  z/L  are  all  less  than  1.  This 
shows  that  the  logarithmically  linear  law  can  be  satisfied  very 
well  up  to  a  10  kilometer  range.  For  type  D,  the  mean  altitude 
of  the  smog  10  kilometers  away  is  190  meters.  We  can  still 
consider  that  it  is  dispersing  in  a  constant  flux  layer  in 
approximation.  However,  for  types  C-A,  the  constant  flux  layer 
is  exceeded. 

Therefore,  for  types  C-A,  due  to  the  limitation  of  the 
constant  flux  layer  hypothesis,  the  similarity  theory  on  the 
ground  level  is  no  longer  satisfied  when  the  dispersion  model  (6) 
is  applied  to  a  10  kilometer  range.  Despite  the  fact  that  there 
was  a  lack  of  theoretical  basis  in  deriving  the  vertical  dispersion 
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1 

* 
i 

model  (6) ,  however,  as  an  empirical  treatment,  it  does  not  seem 
to  cause  any  large  error.  For  example,  the  results  of  types  A-C 
are  in  considerable  agreement  with  the  results  obtained  by  Briggs. 

Therefore,  as  a  whole  we  cam  recommend  the  vertical  dispersion 
model  (6) .  However,  one  must  remember  that  there  is  no  theoretical 
basis  for  types  A-C  * 

IV.  EFFECTIVE  ROUGHNESS  LENGTH  AND  ITS  COMPUTATIONAL  METHOD  /371 

Horizontal  homogeneity  is  an  important  premise  in  the 
similarity  theory.  However,  the  actual  terrain  is  very  complex. 

It  is  very  difficult  to  satisfy  this  hypothesis.  For  an 

r  19) 

inhomogeneous  topography,  the  meaning  of  z0  will  change'  . 

Usually  the  parameter  zQ  is  determined  on  a  mast  or  a  tower.  It 
is  a  local  parameter.  Strictly  speaking,  it  can  only  represent 
the  value  at  the  measuring  point. 

With  regard  to  an  inhomogeneous  terrain,  Fiedler  and 
Panofsky^5*  presented  the  concept  of  "effective  roughness  length." 

It  could  be  derived  from  the  neutral  wind  profile  theory  under  a 
homogeneous  condition  by  using  the  measured  mean  stress  near  the 
ground  in  the  region  under  consideration. 

The  smog  dispersion  test  data  to  be  used  later  in  this 

(18  20 ) 

paper  had  been  obtained  under  various  topographical  conditions  ' 

Therefore,  the  terrains  were  very  rough  and  inhomogeneous.  We 
adopted  the  following  method  to  derive  the  effective  roughness 
length  zQ  for  each  experimental  point  within  a  range  under 
consideration . 

From  (2)  we  can  obtain 

?- 

where  zm  and  zj  represent  the  mean  vertical  positions  at  distances 
of  Xm  and  Xj  downwind  from  the  same  smog  path,  respectively.  zQ 
represents  the  average  roughness  length  of  the  inhomogeneous 
terrain  between  Xg  and  X j .  with  respect  to  each  smog  path ,  \ 

several  values  of  zQ  could  be  calculated.  Then,  the  average  was  | 


i.flai*  -  1  +  J»1 

_  A  U  i-J 

i  Jin  ii  -  1  +  -*] 

>■  *»  IjJ 


(7) 
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obtained.  Finally,  the  values  of  z0  for  many  smog  paths  were 
used  in  order  to  find  the  total  average,  which  is  the  effective 
roughness  length  taken  into  consideration  by  us  in  each  region. 

Based  on  the  aforementioned  method,  the  values  of  zQ  were 
calculated  for  seven  different  terrains,  and  the  results  are 
shown  in  Table  4. 


Table  4.  Calculated  Results  of  Effective  Roughness  Length  ^(m) 

1.  Terrain;  2.  Number;  3.  Number  of  Smog  Paths; 

4.  Calculated  Using  Equation  (7);  5.  Calculated 
Using  Equation  (8);  6.  Local  Roughness 


1.  *  * 

J.  *  # 

j.  *  a  n 

r.  <•)*«■* 

(.  *««a* 

A 

21 

10 

0.4 

0.33 

0.37 

B 

41 

to 

0.1 

0.M 

0.33 

v  C 

40 

4 

0.2 

0.22 

0.(3 

0 

17 

J 

0.2 

— 

— 

E 

27 

10 

0.1 

— 

_ 

F 

13 

3 

0.2 

— 

— 

G 

11 

1 

0.7 

— 

— 

From  these  results  one  can  see  that  the  effective  roughness 
length  is  not  in  agreement  with  the  locally  measured  roughness 
length  (derived  from  the  average  wind  speed  profiles  measured 
on  towers) . 

In  order  to  further  demonstrate  the  representation  and 

reliability  of  the  z_  calculated  based  on  equation  (7) ,  we 

u  (5) 

changed  the  turbulent  form  used  by  Fiedler  and  Panofsky  into 

the  following,  which  we  could  observe  and  measure: 

£«-*«] P  -{u*  X  (l  - +  Zi^]}  l8) 

where  F  is  a  parameter.  Equation  (8)  was  used  to  calculate  the 
smog  dispersion  test  data  in  three  terrains  identical  to  those 
calculated  using  equation  (7) .  The  results  are  also  shown  in 
Table  4.  From  the  table  one  can  see  that  these  two  computational 
methods  are  in  good  agreement  with  each  other. 


V 
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V.  VERTICAL  DISPERSION  MODEL  IN  HOMOGENEOUS  ROUGH  TOPOGRAPHY 
A  10  kilometer  range  vertical  dispersion  model  has  already 
been  derived  for  a  flat  and  homogeneous  topography,  in 
this  section,  we  will  attempt  to  extend  the  aforementioned 
dispersion  model  to  complex  terrains  using  the  effective  roughness 
length  concept  mentioned  in  Section  (IV) .  The  method  is  to 
primarily  replace  zQ  in  (2) -(5)  by  the  effective  roughness  length 
zc.  This  also  means  that  an  inhomogeneous  rough  terrain  is 
considered  to  be  a  homogeneous  terrain  with  an  effective  roughness 
length  zQ  in  the  treatment. 

In  addition,  as  we  pointed  out  before,  m  is  related  to  the 
roughness.  The  rougher  the  ground  is,  the  larger  m  becomes.  In 
order  to  derive  the  relation  between  m  and  Zq,  we  used  a  fact 
which  has  been  proven  in  a  macroscopic  observation,  i.e.  the 
effect  of  an  abrupt  change  in  roughness.  It  decreases  as  the 
distance  increases.  At  a  very  far  distance,  the  effect  of  the 
terrain  can  be  neglected.  Hence,  it  is  possible  to  use  equation 
(5)  to  derive  the  values  of  m  corresponding  to  each  stability 
type  for  various  values  of  zQ.  The  values  of  m  when  zQ  *  1  meter 
are  shown  in  Table  5. 

Table  5.  Values  of  m  when  zQ  *  1  meter. 

1.  Stability  Type 


A  comparison  of  the  values  of  m  in  Table  1  to  those  in 
Table  5  shows  that  as  the  terrain  gets  rougher,  the  absolute 
value  of  m  becomes  significantly  larger. 

In  order  to  provide  the  vertical  dispersion  model  for  rough 
terrains,  we  used  the  z0  *  1  meter  terrain  to  represent  a  rough 
terrain.  Using  Table  5  and  the  values  of  L  in  Table  1,  the 


following  vertical  dispersion  models,  similar  to  those  homogeneous 
rough  terrains  in  (6) ,  could  be  obtained 

/•  a£:  X  —  2.11  X  10'  ■+•  4.05  X  1 0*5*  ^  —  6.16  X 

B  18:  X  -  1.57  X  10'  4-  4.9  X  lO'i**"  —  6.47  X  10'i*-* 

j.  C  35:  X  —  1.12  X  10'  +  7.65  X  10'i**0  —  8.77  X  io‘i,,n 

*.  DS5:  X  -  6.76  +  1.5S  X  10'i'-"  -  1.65  X  10,i'  *B 

S.  E  95:  X  -  7.69ilni  -  8i  +  8.04  +  7.08  X  10-*fJ  +  (9) 

1.54  X  KT'i'Ini  -4-  4.1  X  lO-'i1 
F  38:  X  —  7.69ilni  —  9i  +  9.19  4-  2.17  X  + 

6.59  X  10~*i^ni  +  7.54  x  lO-'i* 


1.  Type  A;  2.  Type  B;  3.  Type  C;  4.  Type  D;  5.  Type  E;  6.  Type  F. 

In  order  to  quantitatively  determine  the  magnitude  of  the 

effect  of  ground  roughness  on  j  ,  we  carried  out  a  comparative 

z 

analysis  on  the  values  of  i  calculated  using  (9)  and  (6) .  Within 

Z 

the  0.1  -  10  kilometer  range,  the  ratios  of  o,  for  the  same 
distance  and  same  stability  type  are  given  in  Table  6.  From  the 
table  one  can  see  that:  (1)  the  effect  of  terrain  roughening 
caused  j  to  increase,  and  this  effect  decreases  with  increasing 

Z 

X  (rising  z) ;  (2)  as  it  varies  from  instability  toward  stability, 
the  effect  of  roughness  also  increases.  Under  a  stable  situation, 
the  effect  of  terrain  is  the  largest,  the  range  affected  is  the 
farthest  and  thickest,  and  the  magnitude  of  increase  for  a  is 

Z 

also  the  largest;  (3)  in  a  neutral  condition,  the  effect  of  a 

rough  terrain  could  cause  az  to  increase  to  2-3  times  of  that 

in  a  flat  terrain  within  1  kilometer.  This  conclusion  is  in  good 

agreement  with  an  actually  measured  dispersion  parameter  (at  800 

( 18 ) 

meters)  in  a  certain  complex  terrain  which  showed  that  on  the 

average  it  was  1.9-2. 5  times  that  in  a  plain. 

In  order  to  verify  the  applicability  of  the  aforementioned 

results  under  a  complex  topography  condition,  we  calculated  the 

vertical  dispersion  parameters  for  seven  complex  terrains, 

(18  20 ) 

respectively  ’  .  Moreover,  a  comparison  was  made  with  respect 

to  the  actual  smog  data.  The  result  showed  that:  for  data 
collected  on  208  occasions  under  neutral  conditions, 


a  measured/  oz  calculated  =  0.9*1  ±  0.22.  Among  108  pieces  of 
information  in  non-neutral  conditions,  cz  measured/  oz  calculated 
-  o.96  ±  0.21.  Such  results  are  considered  to  be  satisfactory. 

Table  6.  Ratios  of  az  1  for  Z Q  =  1  meter  to  a2  Q<cl  for  ZQ  =  /373 

0.01  Meter 

1.  Stability  Type;  2.  X  (meter) . 


'Vsvaa,,  i 

A 

B 

— 

c 

1 

0 

E 

F 

;oo 

1.41 

1.74 

2.14 

7.81 

2.98 

3.18 

150 

1.24 

1.55 

1.96 

2. 48 

275 

2.88 

200 

1.16 

1.4< 

1.12 

2.31 

2.61 

2  51 

500 

1.02 

1.2i 

1.48 

1.95 

2.08 

1.96 

700 

0.9J 

1.16 

1.41 

1.78 

1.97 

1.79 

1000 

0  95 

1.13 

1.36 

1.6  6 

1.85 

i  174  • 

2000 

0.91 

1.05 

1.24 

1.45 

1.66 

1.59 

1 

5000 

,  0.52 

1  1.00 

1.10 

1.26 

. 

1.44 

| 

7000 

0.92 

1  1.00 

1.08 

1.20 

1 .44 

1-33 

104 

0  92 

1.00 

1.06 

1.15 

1.39 

1.15 

The  experimental  data  used  in  the  earlier  verification  were 
obtained  with  a  sampling  distance  ranging  from  40-520  meters. 

Due  to  the  limitation  of  the  photographic  method,  the  sampling 
distance  is  too  short.  In  order  to  analyze  the  applicability  of 
the  results  obtained  in  this  work  at  a  farther  distance,  we 
carried  out  a  comparative  analysis  using  the  dispersion  test  data 

( 21) 

obtained  in  the  city  of  St.  Louis  in  the  United  States  of  America 
(test  data  on  12  occasions  in  the  evening  under  neutral  conditions 
at  a  sampling  distance  ranging  from  0.78-7.8  kilometers) . 

Pasquill^21^  used  the  following  formula  to  match  some  of  the 
test  data: 


2 


+  -j*  (1  — *  Inc)  —  1  j 


(10) 


-10- 


where  c  *  0.6,  zQ  »  3  meters.  We  used  equations  (2)  and  (3) 
separately  to  conduct  the  corresponding  calculation.  The  ratios 
of  the  calculated  values  to  the  actual  values  are  shown  in  Table  7. 
For  the  purpose  of  comparison,  the  calculated  results  obtained  by 
Pasquill  are  also  given  in  the  table. 

Table  7 .  Comparison  of  Calculated  and  Actually  Measured 
St.  Louis  Data  under  Neutral  Conditions. 

1.  Calculation  Formula;  2.  (2)  z  =  3  meters; 

3.  (3)  z0  *  3  meters  -  L  =  104  miters; 

4.  (3)  z_  =  3  meters  -  L  =-104  meters; 

5  1  /t 

calculated  measured 


From  these  results  one  can  see  that  the  three  models  adopted 
by  us  are  more  consistent  with  the  reality  than  (10)  within  the 
8  kilometer  range. 

In  order  to  further  analyze  the  reliability  of  (9) ,  we  carried 

out  a  comparative  analysis  on  our  results  and  the  corresponding 

(22) 

results  obtained  by  Smith  (derived  from  the  two-dimensional 

dispersion  equation  numerically) .  The  ratios  of  a  within  the 

2 

range  from  0.1-10  kilometers  are  given  in  Table  8. 

Table  8.  Comparison  of  Results  Obtained  by  (9)  and  by  Smith  , 
in  0.1  -  10  Kilometer  Range. 

1.  Stability  Type;  2.  a  (this  work) /a  (Smith) 

Z  2 


From  the  table  one  can  see  that  the  two  are  in  good 
agreement  with  each  other.  With  the  exception  of  Type  B,  the 
ratio  between  the  two  is  always  within  a  factor  of  2.  The 
consistency  is  the  best  under  neutral  conditions.  Moreover,  the 


scatter  is  also  very  small.  But  for  types  C-A,  the  results 
obtained  in  this  work  are  larger  than  those  obtained  by  Smith. 

For  types  E  and  F,  the  results  obtained  by  this  work  are  slightly 
larger.  In  summary,  although  there  are  some  differences  between 
the  two  vertical  dispersion  models  for  homogeneous  rough  terrain 
derived  by  two  different  methods,  the  consistency  of  the  two 
is  still  very  good  from  the  practical  point  of  view. 

VI.  SUMMARY  AND  CONCLUSIONS 

Through  the  analysis  above,  it  is  possible  to  obtain  the 
following  viewpoints: 

1.  For  a  flat  and  homogeneous  terrain,  the  effectiveness 
of  using  an  empirical  method  based  on  the  similarity  theory 
concept  to  estimate  the  vertical  dispersion  within  a  10  kilometer 
range  is  equivalent  to  that  of  the  Briggs  intrapolation  formula 
which  is  widely  used  and  commonly  recognized  at  the  present 
moment.  This  dispersion  model  could  be  conveniently  extended  from 
a  homogeneous  flat  terrain  to  a  homogeneous  rough  terrain. 

However,  in  an  unstable  situation,  due  to  the  limitation  of  the 
constant  flux  layer,  there  is  still  a  lack  of  theoretical  basis 
despite  the  fact  that  it  is  possible  to  extend  it  to  a  10 
kilometer  range  based  on  experience.  This  would  require  the 
further  development  of  a  planetary  boundary  layer  similarity 
theory . 

2.  For  a  complex  terrain,  the  similarity  theory  could  not 
be  used,  in  principle.  However,  in  reality,  it  is  not  possible 
to  really  find  a  completely  homogeneous  terrain.  Essentially, 
it  is  necessary  to  consider  -the  extent  of  the  effect  of  this 
inhomogeneity  on  the  application  of  the  theory.  In  this  paper, 
the  relation  between  X  and  z  in  a  neutral  condition  was  used  to 
determine  the  effective  roughness  length.  The  results  were 
compared  with  those  obtained  using  the  method  of  Fiedler  and 
Panofsky,  which  showed  that  the  method  used  in  this  work  was 
feasible. 

Using  the  effective  roughness  length  obtained  through 
calculation  in  this  work,  the  dispersion  model  of  a  homogeneous 
flat  terrain  was  extended  to  rough  terrains.  A  comparison  with 


the  actually  observed  data  and  the  relevant  results  obtained  by 
Smith  showed  a  very  good  agreement. 

3.  A  comparison  of  the  vertical  dispersion  characteristics 
of  various  rough  terrains  showed  that  the  rougher  the  terrain  is, 
the  larger  the  dispersion,  parameter  becomes.  For  the  same 
roughness,  the  effect  of  terrain  on  dispersion  gradually  increases 
from  instability  to  stability.  It  is  most  appaxent  in  stable 
situations.  The  actual  measured  data  also  supported  this  fact'  '  ' 
In  summary,  for  homogeneous  flat  terrains,  from  a  neutral 
to  a  stable  layer,  it  is  possible  to  use  the  similarity  theory 
concept  to  empirically  estimate  the  vertical  dispersion  in  10 
kilometers.  For  complex  terrains,  after  replacing  local  roughness 
lengths  with  effective  roughness  lengths,  the  similarity  theory 
can  also  be  applied  in  approximation.  From  this,  an  empirical 
vertical  dispersion  model  could  be  obtained/  With  regard  to  its 
applications  in  the  areas  of  atmospheric  pollution,  environmental 
quality  evaluation,  and  engineering  design,  the  accuracy  is 
sufficient. 
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